Quark fragmentation to n"^, 7r°, K"^, p and p in the nuclear environment 

A. Airapetian,'^'' N. Akopov,^° Z. Akopov,'^" M. Amarian,^' ^° V.V. Ammosov,^^ A. Andrus/^ E.G. Aschenauer,^ 
W. Augustyniak,^^ R. Avakian,^° A. Avetissian,^° E. Avetissian/° P. Bailey/^ V. Baturin,-^^ C. Baumgarten,^^ 
M. Beckmann,^ S. Belostotski,^! S. Bernreuther,^ N. Bianchi/° H.P. Blok,20'28 H. Bottcher,^ A. Borissov," 
A. Borysenko,!" M. Bouwhuis,!^ J. Brack,* A. BriiU/^ V. Bryzgalov,22 G.P. Capitani,i" B.C. Chiang,!^ G. Ciullo,^ 
M. Gontalbrigo,^ P.F. Dalpiaz,^ R. De Leo,^ L. De Nardo,i E. De Sanctis/° E. Devitsin,!^ P. Di Nezza/o 
M. Duren/3 Ehrenfried,^ A. Elalaoui-Moulay,^ G. Elbakian,30 F. EUinghaus,^ U. Elschenbroich,!^ j ^iy4 
R. Fabbri,^ A. Fantoni,io A. Fechtchenko/ L. Felawka,^^ B. Fox,"* J. Franz," S. Frullani,^* G. Gapienko,^^ 
V. Gapienko,22 F. Garibaldi,^* K. Garrow,i^25 Garutti,20 D. Gaskell,* G. Gavrilov,5.26 y. Gharibyan,^" 
G. Graw,i9 O. Grebeniouk,^! L.G. Greeniaus.^-^^ I.M. Gregor,^ K. Hafidi,^ M. Hartig,^^ D. Hasch,!*' 
D. Heesbeen,20 M. Henoch,^ R. Hertenberger,!^ W.H.A. Hesselink,20'28 A. Hillenbrand,^ M. Hoek,i3 Y. Holler,^ 
B. Hommez,^^ G. larygin,'^ A. Ivanilov,^^ A. Izotov,^^ H.E. Jackson,^ A. Jgoun,^^ R. Kaiser,^"* E. Kinney,* 
A. Kisselev,^^ K. Konigsmann," M. Kopytin,^ V. Korotkov,^ V. Kozlov,^^ B. Krauss,® V.G. Krivokhijine,^ 
L. Lagamba,'^ L. Lapikas,^° A. Laziev,^°'^* P. Lenisa,^ P. Liebing,^ T. Lindemann,^ K. Lipka,^ W. Lorenzon,^'' 
J. Lu,26 B. Mailieu,i2 N.G.R. Makins,!^ B. Marianski,^^ H. Marukyan,^'' F. Masoli,^ V. Mexner,^^ 
N. Meyners,^ O. Mikloukho,^^ C.A. Miller, ^'^^ Y. Miyachi,^^ V. Muccifora,^" A. Nagaitsev,^ E. Nappi,^ 
Y. Naryshkin,2i A. Nass,** M. Negodaev,^ W.-D. Nowak,*^ K. Oganessyan,^' H. Ohsuga,^^ N. Pickert,^ 
S. Potasliov,!^ D.H. Potterveld,^ M. Raithel,** D. Reggiani,^ P.E. Reimer,^ A. Reischl,20 A.R. Reolon,i° 

C. Riedl,^ K. Ritli,** G. Rosner,** A. Rostomyan,30 L. Rubacek," D. Ryckbosch,!^ Y. Salomatin,22 
I. Sanjiev,2.2i l. Savin,^ C. Scarlett," A. Schafer,23 G. SchiU," G. Schnell,^ K.P. Schiiler,^ A. Schwind,^ 
J. Seele,^^ R. Seidl,** B. Seitz,!^ R. Slianidze,** G. Shearer, i* T.-A. Shibata,^^ V. Shutov,'' M.C. Simani,20. 28 
K. Sinram,^ M. Stancari,^ M. Statera,^ E. StefFens,^ J.J.M. Steijger,^^ H. Stenzel,^^ J. Stewart,^ U. Stosslein,* 
P. Tait,* H. Tanaka,27 S. Taroian,30 B. Tcliuiko,^^ A. Terkulov,!^ A. Tkabladze,!^ a. Trzcinski,^^ 
M. Tytgat,i2 a. Vandenbroucke,!^ ^^j^ (jgj. j^f^^^ 20,28 q ^^j. Steenhoven,^^ M.G. Vetterli,^^' 26 
V. Vikhrov,2i M.G. Vincter,i G. Vogel,** M. Vogt.^ J. Volnier,^ C. Weiskopf,** J. Wendland,25. 26 
J. Wilbert,^ G. Ybeles Smit,28 S. Yen,26 B. Zihlmann,20' 28 Zolirabian,^^ and P. Zupranski29 

(The HERMES GoUaboration) 

^Department of Physics, University of Alberta, Edmonton, Alberta T6G 2J1, Canada 
'"^Physics Division, Argonne National Laboratory, Argonne, Illinois 60439-4843, USA 
Istituto Nazionale di Fisica Nucleare, Sezione di Bari, 70124 Bari, Italy 
^Nuclear Physics Laboratory, University of Colorado, Boulder, Colorado 80309-0446, USA 
'^DESY, Deutsches Elektronen- Synchrotron, 22603 Hamburg, Germany 
^DESY Zeuthen, 15738 Zeuthen, Germany 
'Joint Institute for Nuclear Research, 141980 Dubna, Russia 
^ Physikalisches Institut, Universitdt Erlangen-Niimberg, 91058 Erlangen, Germany 
^Istituto Nazionale di Fisica Nucleare, Sezione di Ferrara and 
Dipartimento di Fisica, Universita di Ferrara, 44^00 Ferrara, Italy 
^"istituto Nazionale di Fisica Nucleare, Laboratori Nazionali di Frascati, 00044 Frascati, Italy 

"Fakultat fiir Physik, Universitdt Freiburg, 79104 Freiburg, Germany 
^''^ Department of Subatomic and Radiation Physics, University of Gent, 9000 Gent, Belgium 
''^ Physikalisches Institut, Universitdt Giefien, 35392 GiefSen, Germany 
^'^ Department of Physics and Astronomy, University of Glasgow, Glasgow G12 8QQ, United Kingdom 
Department of Physics, University of Illinois, Urbana, Illinois 61801-3080, USA 
^''Laboratory for Nuclear Science, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA 
' 'Randall Laboratory of Physics, University of Michigan, Ann Arbor, Michigan 48109-1120, USA 
Lebedev Physical Institute, 117924 Moscow, Russia 
'^Sektion Physik, Universitdt Miinchen, 85748 Carching, Germany 
Nationaal Instituut voor Kernfysica en Hoge-Energiefysica (NIKHEF), 1009 DB Amsterdam, The Netherlands 
'^'^ Petersburg Nuclear Physics Institute, St. Petersburg, Gatchina, 188350 Russia 
'^''Institute for High Energy Physics, Protvino, Moscow region, 142281 Russia 
^^Institut fiir Theoretische Physik, Universitdt Regensburg, 93040 Regensburg, Germany 
^'^Istituto Nazionale di Fisica Nucleare, Sezione Roma 1, Gruppo Sanitd 
and Physics Laboratory, Istituto Superiore di Sanitd, 00161 Roma, Italy 
'^^ Department of Physics, Simon Eraser University, Burnaby, British Columbia V5A 1S6, Canada 
'^"TRIUMF, Vancouver, British Columbia V6T 2A3, Canada 
''^''Department of Physics, Tokyo Institute of Technology, Tokyo 152, Japan 



2 



Department of Physics and Astronomy, Vrije Universiteit, 1081 HV Amsterdam, The Netherlands 
'^^ Andrzej Soltan Institute for Nuclear Studies, 00-689 Warsaw, Poland 
'^''Yerevan Physics Institute, 375036 Yerevan, Armenia 

The influence of the nuclear medium on lepto-production of hadrons was studied in the HERMES 
experiment at DESY in semi-inclusive deep-inelastic scattering of 27.6 GeV positrons off deuterium, 
nitrogen and krypton targets. The differential multiplicity for krypton relative to that of deuterium 
has been measured for the first time for various identified hadrons (tt^, tt", tt", K'^ , K~ , p and p) 
as a function of the virtual photon energy u, the fraction z of this energy transferred to the hadron, 
and the hadron transverse momentum squared pf. The multiplicity ratio is strongly reduced in 
the nuclear medium at low i/ and high z, with significant differences among the various hadrons. 
The distribution of the hadron transverse momentum is broadened towards high pt in the nuclear 
medium, in a manner resembling the Cronin effect previously observed in collisions of heavy ions 
and protons with nuclei. 



PACS numbers: 13.87.Fh, 13.60.-r, 14.20.-c, 14.40.-n gluon-bremsstrahlung model for leading hadron atten 



The phenomenon of confinement in QCD imposes it- 
self dynamically in the process of hadronization, i.e. the 
mechanism by which final-state hadrons are formed from 
a quark that has been struck hard. This process, also 
known as the fragmentation of quarks into hadrons, can 
be described by fragmentation functions Dj{z), denot- 
ing the probability that a quark of flavor / produces a 
hadron of type h carrying a fraction z of the energy of 
the struck quark in the target rest frame. In the nuclear 
medium additional soft processes may occur before the 
final-state hadron is completely formed. The nuclear en- 
vironment may thereby influence the hadronization pro- 
cess, e.g. cause a change in the quark fragmentation 
functions, in analogy to the EMC finding of a medium 
modification of the quark distribution functions 0|. 

The understanding of quark propagation in the nu- 
clear medium is crucial for the interpretation of ultra- 
relativistic heavy ion collisions, as well as high en- 
ergy proton- nucleus and lepton- nucleus interactions 0. 
Quark propagation in the nuclear environment involves 
processes like multiple interactions with the surrounding 
medium and induced gluon radiation, resulting in energy 
loss of the quark. If the final hadron is formed inside the 
nucleus, the hadron can interact via the relevant hadronic 
interaction cross section, causing a further reduction of 
the hadron yield. Therefore, quark and hadron propaga- 
tion in nuclei are expected to result in a modification i.e. 
a "softening" of the leading-hadron spectra Q compared 
to that from a free nucleon. By studying the properties 
of the leading-hadrons emerging from nuclei, information 
on the characteristic time-distance scales of hadroniza- 
tion can be derived. 

The hadronization process in the nuclear medium is 
traditionally described in the framework of phenomeno- 
logical string models 0, 0, ^^'^ final state inter- 
actions of the produced hadrons with the surrounding 
medium |^. Alternatively, in- medium modifications of 
the quark fragmentation functions have been proposed, 
either expressed in terms of their nuclear rescaling 
or parton energy loss [U and hi gher -twist contributions 
to the fragmentation functions [iJI, or in terms of a 
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uation [l^j. The models of Refs. pH also incorpo- 
rate hadronic final state interactions. These recent QCD- 
inspired models provide a theoretical description of the 
hadronization process in deeg-inelastic scattering, rela- 
tivistic heavy-ion collisions and Drell-Yan reactions 
on nuclear targets [llLll^. Moreover, some of these mod- 
els contain a so-far untested QCD prediction, i.e. that 
the induced radiative energy loss of a quark traversing a 
length L of hot or cold matter is proportional to Il5l 
due to the coherence of the gluon radiation process jl6j . 

Semi-inclusive deep-inelastic lepton-nucleus collisions 
are most suitable to obtain quantitative information on 
the hadronization process. In contrast to hadron-nucleus 
and nucleus-nucleus scattering, in deep-inelastic scatter- 
ing no deconvolution of the parton distributions of the 
projectile and target particles is needed, so that hadron 
distributions and multiplicities from various nuclei can 
be directly related to nuclear effects in quark propaga- 
tion and hadronization. 

The experimental results for semi-inclusive deep- 
inelastic scattering on nuclei are usually presented in 
terms of the hadron multiplicity ratio -Rj^, which rep- 
resents the ratio of the number of hadrons of type h pro- 
duced per deep-inelastic scattering event on a nuclear tar- 
get of mass A to that from a deuterium target (D). The 
ratio R'^i depends on the leptonic variables v and Q^, 
where v and —Q^ are the energy in the target rest frame 
and the squared four- momentum of the virtual photon re- 
spectively, and on the hadronic variables z and , where 
Pt is the hadron momentum component transverse to the 
virtual photon direction. Fig. ^ illustrates the definition 
of the relevant lepton and hadron kinematic variables for 
this analysis. The multiplicity ratio is defined as: 



RIi{z,p,pIQ^ 



(1) 



where is the yield of semi-inclusive hadrons in a 
given (z, i^,Pf,Q^)-hm, and the yield of inclusive deep- 
inelastic scattering leptons in the same (i/,Q^)-bin. The 
ratio R'Ij is usually evaluated as a function of and z 
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FIG. 1: Kinematic planes for hadron production in semi- 
inclusive deep-inelastic scattering and definitions of the rel- 
evant lepton and hadron variables. The quantities k {k') and 
E {E') are the 4-momentum and the energy of the incident 
(scattered) positron, and ph and Eh are the 4-momentum and 
the energy of the produced hadron. 

only, while integrating over all other kinematic variables, 
as existing data for R^^ show a weak dependence on ei- 
ther Q2 orp2 

In the past, semi-inclusive leptoproduction of undif- 
ferentiated hadrons from nuclei was studied at SLAC 
with electrons JJJ, and at CERN and FNAL with high- 
energy muons by EMC and E665 E(Tj | . Recently, 
HERMES reported more precise data [2l| on the pro- 
duction of charged hadrons, as well as identified tt"*" and 
7r~ mesons, in deep-inelastic positron scattering on ni- 
trogen relative to deuterium. A significant difference was 
found between the multiplicity ratio of positive and neg- 
ative hadrons, while the multiplicity ratio for identified 
pions was found to be the same for both charges. In or- 
der to clarify this issue, additional measurements of R'^j 
with identification of other hadron species have been per- 
formed at HERMES. 

In this paper we present results on the hadron multi- 
plicities on krypton relative to deuterium, providing the 
first measurements of the multiplicity ratio for identified 
pions, kaons, protons and antiprotons. Additionally the 
nitrogen data for charged hadrons and identified pions 
are reevaluated, now covering a wider kinematic range 
than in Ref. \2l\ . The measurements described were per- 
formed with the HERMES spectrometer using the 
27.6 GeV positron beam stored in the HERA ring at 
DESY. The spectrometer consists of two identical halves 
located above and below the positron beam pipe. Both 
the scattered positrons and the produced hadrons were 
detected simultaneously within an angular acceptance of 
± 170 mrad horizontally, and ± (40 - 140) mrad verti- 
cally. 

The data were collected using krypton and deuterium 
gas targets internal to the positron storage ring. Either 
polarised deuterium or unpolarised high density kryp- 
ton gas was injected into a 40 cm long tubular open- 



ended storage cell. Target areal densities up to 1.4 x 
10^^ nucleons/cm^ were obtained for krypton. During 
these high-density runs HERA operated in a dedicated 
mode for the HERMES experiment. This made it pos- 
sible to accumulate the krypton statistics within a few 
days in 1999. The deuterium data were collected over a 
period of one year (1999) with a lower-density polarised 
target. The yields from deuterium were averaged over 
the two spin orientations. 

The positron trigger was formed by a coincidence be- 
tween signals from three scintillator hodoscope planes, 
and a lead-glass calorimeter where a minimum energy 
deposit of 3.5 GeV (1.4 GeV) for unpolarised (polarised) 
target runs was required. The identification of the scat- 
tered positrons was accomplished using a transition- 
radiation detector, a scintillator preshower counter, and 
an electromagnetic calorimeter. 

The identification of charged pions, kaons, protons 
and antiprotons was accomplished using the information 
from the RICH detector |23| , which replaced a threshold 
Cerenkov counter used in the previously reported mea- 
surements on ^^N [i^. This detector uses two radia- 
tors, a 5 cm thick wall of silica areogel tiles followed by a 
large volume of C4F10 gas, to provide separation of pions, 
kaons, and protons over most of the kinematic acceptance 
of the spectrometer. The pions and kaons identified by 
the RICH detector are analysed in the momentum re- 
gion between 2.5 GeV and 15 GeV, while for the iden- 
tified protons and antiprotons the momentum region is 
restricted to the range between 4 GeV and 15 GeV in 
order to reduce possible contaminations from misidenti- 
fied hadrons. The identification efficiencies and contam- 
inations for pions, kaons, protons and antiprotons have 
been determined in a Monte Carlo simulation as a func- 
tion of the hadron momentum and multiplicity in the 
relevant detector half. These RICH performance param- 
eters were verified in a limited kinematical domain using 
known particle species from identified resonance decays. 
They were used in a matrix method to unfold the true 
hadron distributions from the measured ones. 

The electromagnetic calorimeter ji^l provided neutral 
pion identification by the detection of two neutral clus- 
ters originating from the two decay photons. Each of the 
two clusters was required to have an energy E-y > 0.8 
GeV. The background was evaluated in each kinematic 
bin by fitting the two-photon invariant mass spectrum 
with a Gaussian plus a polynomial that fits the shape of 
the background due to uncorrelated photons. The num- 
ber of detected tt'^ mesons was obtained by integrating 
the peak, corrected for background, over the ±2ct range 
with respect to the centroid of the Gaussian. The identi- 
fied neutral pions were analysed in the same momentum 
range as the charged pions, i.e. between 2.5 GeV and 15 
GeV. 

Scattered positrons were selected by imposing the con- 
straints Q2 > 1 GeV^, W = .j2Mv + A'P > 2 GeV 
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for the invariant mass of the photon-nucleon system 
where M is the nucleon mass, and y = v/E < 0.85 
for the energy fraction of the virtual photon. The re- 
quirements on W and y are applied to exclude nucleon 
resonances and to limit the magnitude of the radiative 
corrections, respectively. In the previously reported data 
on ^^N [i^l an additional constraint x = > 0.06 

for the Bjorken scaling variable was applied, in order to 
exclude the kinematic region where an anomalous ratio 
of the longitudinal to transverse cross sections for deep- 
inelastic scattering from ^^N appeared. Since that report 
this anomaly was found to be due to a peculiar local in- 
strumental inefficiency, for which corrections have now 
been evaluated [i^. By applying these corrections the 
cc-range was extended down to x=0.02 both for the new 
data and for the previously published ^'^N data. 

Most of the hadrons from target fragmentation were 
suppressed by the requirement z > 0.2. The lower 
constraint on the hadron momentum implies that the 
2;— acceptance is restricted to rather large values as v de- 
creases. Hence, to ensure that the z^-dependence of the 
multiplicity ratio does not correspond to a strong varia- 
tion of the mean z-value, the present data were confined 
tov>7 GeV. 

Under the kinematic constraints described above, the 
number of selected deep-inelastic scattering events is 8.2 
(7.3) xlO^ for krypton (deuterium). The number of 
charged pions is 11.3 (13.5) xlO*, while the number of 
neutral pions is 1.9 (2.3) xlO''. The number of kaons is 
2.1 (2.3) xlO**, and the number of protons and antipro- 
tons is 1.1 (1.1) xlO''. 

The data have been corrected for radiative processes 
involving nuclear elastic, quasi-elastic and inelastic scat- 
terin g, u sing the codes of Refs. [2^ The code of 

Ref. [231 was modified to include the measured semi- 
inclusive deep-inelastic scattering cross sections. The size 
of the radiative corrections applied to R'^^ was found to 
be negligible in most of the kinematic range, with a max- 
imum of about 7% at the highest value of v, as most of 
the radiative contributions cancel in the multiplicity ra- 
tio 

The charged pion sample is contaminated by pions 
originating from the decay of heavier mesons. The main 
effect on the pion multiplicities is due to the decay of 
exclusively produced vector mesons, which may affect 
the multiplicities by an amount ranging from about 1% 
at low z up to 30% (45%) at high z for positive (nega- 
tive) pions, as has been estimated from a Monte Carlo 
simulation. The effect on the super-ratio R'Ij is smaller, 
but does not cancel as the p° vector meson are also at- 
tenuated in the nuclear medium. By taking into account 
the measured p*^ nuclear transparency 29] , the remaining 
effect on R'Ij has been estimated and included in the sys- 
tematic uncertainty. No correction was applied to R'l^. 
Pions resulting from the decay of p° mesons formed in 
the fragmentation process are included in R'Ij- 



The systematic uncertainty is reduced due to the fact 
that super-ratios of seminclusive and inclusive yields are 
measured. The contributions to the systematic uncer- 
tainty of R'Ij arise from radiative corrections (< 2%), 
hadron identification (1.5% for neutral pions, 0.5% for 
kaons, 1% for protons and 2% for antiprotons), overall 
efficiency (< 2%), and p°-meson production for positive 
(0.3% - 4%) and negative (0.3% - 7%) pions 

The geometric acceptance for semi-inclusive hadron 
production has been verified to be the same for both the 
krypton and deuterium targets by studying the multiplic- 
ity ratio as a function of the hadron polar angle. This 
ratio was found to be constant within the experimental 
precision. 

The multiplicity ratio has been determined as a func- 
tion of either z, v 01 p^, while integrating over all other 
kinematic variables. In Fig.|21the multiplicity ratios for 
all charged hadrons with z > 0.2 are presented as a func- 
tion of i> together with data of previous experiments on 
nuclei of similar size. In the top panel the HERMES data 
on Kr are compared with the SLAC [3l and CERN ^ 
data for Cu. In the lower panel the reevaluated HER- 
MES data on ^^N are displayed together with data on 
[ll[l3. Due to the extension down to a;=0.02, the 
^'^N data shown in Fig. |21 have a higher statistical accu- 
racy than the data reported in Ref. [13 . The HERMES 
data for i?^ are observed to increase with increasing 
roughly approaching the EMC results at higher values of 
v. The discrepancy with the SLAC data is partially due 
to the fact that semi-inclusive cross section ratios were 
measured at SLAC instead of multiplicity ratios. Thus no 
corrections in the SLAC data were made for the target- 
mass dependence of the inclusive deep- inelastic scattering 
cross section, as discussed in Ref. 21]. 

A stronger attenuation is observed for Kr than for ^^N, 
the average ratios and the total experimental uncertain- 
ties being i?5^=0.802 ± 0.021 and i?^,^=0.954 ± 0.023, re- 
spectively. Fig.|3shows the dependence on z of the same 
multiplicity ratios for v > 7 GeV. This figure includes 
the region z < 0.2, which contains contributions from 
both target fragmentation hadrons and leading hadrons 
decelerated in nuclear re-scattering. Qualitatively, the 
dependences on ly and z of the Kr data resemble those of 
the "'^^N data, but the features are more pronounced. 

The measured v- and z-dependences for both krypton 
and nitrogen are compared to several model calculations 
shown in Figs. E] and El as solid curves 0, and dot- 
dashed curves the nuclear modification 
of hadron production in deep-inelastic scattering is de- 
scribed as a rescaling of the quark fragmentation func- 
tions, supplemented by nuclear absorption. In this model 
the nuclear absorption contribution is dominated by the 
string interaction, while the subsequent interaction of the 
fully formed hadron contributes only a few percent to R% 
for krypton. The calculation overestimates the ^^N at- 
tenuation, but gives a fairly good account of both the 
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FIG. 2: Charged hadron multiplicity ratio R\.[ as a func- 
tion of V for z > 0.2. In the upper panel HERMES data on 
Kr are compared to SLAC |l^ and CERN jl^ data on Cu. 
In the lower panel the HERMES data on are compared 
with CERN and SLAC data on ^^C. The error bars repre- 
sent the statistical uncertainties, and the systematic uncer- 
tainty is shown as the band. The solid curves are calculations 
from Ref. [To| and the dot-dashed curves are calculations from 
Ref. m. 



V and z dependences, for z > 0.2, of the Kr data. In 
Ref. JM nuclear modification of the quark fragmenta- 
tion process in deep-inelastic scattering has been evalu- 
ated taking into account multiple parton scattering and 
induced energy loss in the medium. The only free param- 
eter of this model was tuned to reproduce the ^''N data, 
and the derived energy loss has been used to determine 
the initial gluon density in the Au-f Au collision data col- 
lected by the PHENIX collaboration at RHIC 113 . The 
model roughly reproduces the changes in R'^.j when going 
from nitrogen to krypton. 

For the first time the I'-dependences of the multiplicity 
ratio were studied for identified neutral and charged pi- 
ous, kaons, protons and antiprotons, as shown in the left 
part of Fig. ^ The corresponding z-dependences of R'^yf 
with V > 7 GeV are shown in the right part of Fig. ^ In 
the bottom panels the average values for and z or v 



FIG. 3; Charged hadron multiplicities Rf^j as a function of 
z for u > 7 GeV. The data are compared to the same calcu- 
lations as in Fig. |5| The error bars represent the statistical 
uncertainties and the systematic uncertainty is shown as the 
band. 

are displayed for all the presented data. The results for 
both charge states of the pion and the kaon are compared 
to the calculations of Ref. 0. A good description for 
the TT* and data is observed, while the attenuation 
is under-predicted for the K~ data. No model predic- 
tions are available for protons and antiprotons. Average 
R'li values are obtained by integrating yields over ly and 
z. The results presented in Fig. 01 and the average R'lf 
values reported in TablelUfor z > 0.2, show that the mul- 
tiplicity ratios for positive and negative pions are similar, 
in agreement with what was already found on -'^'^N . In 
addition, the multiplicity ratio for neutral pions is found 
to be consistent, within the total experimental uncertain- 
ties, with that for charged pions as well as for negative 
kaons. However, R'^j for positive kaons is significantly 
larger. An even larger difference is observed between 
protons and their antiparticles compared to the meson 
case. These differences in R'Ij of positive and negative 
kaons, as well as those between protons and antiprotons, 
are still present at higher values of z. This is shown in 
the last column of Tabled where the average R'Ij values 
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FIG. 4: Multiplicity ratios for identified pions, kaons, protons and antiprotons from a Kr target as a function of v for z > 0.2 
(left), and as a function of z for v > 1 GeV (right). In the upper right panel the multiplicity ratio for identified pions from a 
^''N target are also shown. The closed (open) symbols represent the positive (negative) charge states, and the crosses represent 
tt" mesons. In the bottom panels the average z and v values are displayed: pions and kaons (protons and antiprotons) are 
shown as closed (open) circles; the average values are indicated by the open stars referring to the right-hand scales. The 
inner (outer) error bars represent the statistical (total) uncertainties. The thick (thin) solid curves represent the calculations 
of Ref. for positive (negative) charge states. Multiplicity ratios for negative kaons and antiprotons at the highest z-bins 
are not displayed due to their poor statistical significance. 
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are reported for z > 0.5, i.e. when emphasising leading 
hadrons. In addition the z > 0.5 range is most suitable 
to compare R'lf of mesons and baryons as this compari- 
son is performed at the same average v as shown in the 
bottom right panel of Fig. ^ 



1.5 



h-type 


< Rm > 
z > 0.2 


< i^M > 

z > 0.5 


7V + 


0.775 ± 0.019 


0.712 ±0.023 




0.770 ± 0.021 


0.731 ±0.031 


^« 


0.807 ± 0.022 


0.728 ± 0.024 


K+ 


0.880 ±0.019 


0.766 ± 0.024 


K- 


0.783 ±0.021 


0.668 ± 0.036 


P 


0.977 ± 0.027 


0.816 ±0.029 


P 


0.717 ±0.038 


0.705 ± 0.067 



TABLE I: Multiplicity ratio values for krypton and deuterium 
yields integrated over z and f > 7 GeV. Total experimental 
uncertainties are quoted. 



in< 



It has been suggested that the observed differences 
can be attributed to the mixing of quark and gluon 
fragmentation functions. This mixing gives a different 
modification of the quark and antiquark fragmentation 
functions in nuclei, thus leading to a more significant 
difference between the multiplicity ratio of protons and 
antiprotons than between those of mesons. The observed 
differences in the multiplicity ratios can also be inter- 
preted in terms of different formation times of baryons 
and mesons |3l( . or in terms of different hadron-nuclcon 
interaction cross sections [s^- While this cross section 
is similar for positive and negative pions, it is larger for 
negative kaons as compared to positive kaons, and even 
larger for antiprotons than protons, in qualitative agree- 
ment with the trend shown by the data. 

The data for identified charged pions produced on Kr 
and ^"^N for z > 0.5 were used to estimate the mass- 
number dependence of the nuclear attenuation. The ni- 
trogen data, measured in a more restricted momentum 
range between 4-13.5 GeV due to use of the Cerenkov 
detector, are shown in the upper right panel of Fig. 0] 
By using the same constraints on the krypton data and 
assuming a simple v4"-dependcncc of the nuclear attenua- 
tion 1 — Rpj , the experimental data are found to be closer 
to the -dependence predicted in Ref. than 
the A 3 -dependence that follows from models based on 
nuclear absorption effects only. Data on more nuclei are 
needed to enable systematic studies of the ^-dependence 
of the nuclear attenuation effects. 

The pt distribution of the observed hadrons is expected 
to be broadened on a nuclear target compared to a pro- 
ton target due to multiple scattering of the propagating 
quark and hadron. This effect is known as the Cronin 
effect |33j and has previously been observed in heavy- 
ion and hadron- nucleus induced reactions. A nuclear en- 
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FIG. 5: Multiplicity ratio for charged hadrons versus pt for 
f >7 GeV and z >0.2. The HERMES data on Kr and are 
compared to the EMC ^3 data for Cu in the range 10 < ly < 
80 GeV. The error bars represent the statistical uncertainties. 
The systematic uncertainty for Kr (^''N) is shown as the lower 
(upper) band. 



hancement at high is also observed in the HERMES 
data shown in Fig. El both for ^^N and Kr nuclei. In 
this plot the EMC 18] data on Cu which cover a differ- 
ent ly-range, 10 < < 80 GeV, are also displayed. The 
data for < 0.7 GeV^ show the attenuation previously 
discussed, while the data for p^ > 0.7 GeV^ refiect the pt- 
broadening ascribed to multiple scattering effects. This 
effect is similar to the one reported for proton-nucleus 
and nucleus-nucleus collisions ;34 | but is smaller in mag- 
nitude. The enhancement is also predicted to occur at 
a pt-scale of about 1-2 GeV jlM H^l i in agreement with 
the semi-inclusive deep-inelastic scattering data shown 
in Fig. El The HERMES data may help to interpret 
the new relativistic heavy-ion results from SPS "37^ and 
RHIC j3Q], which show a weaker pt enhancement than 
expected from the original Cronin effect. 

In summary, the multiplicities of charged hadrons and 
of identified pions, kaons, protons and antiprotons on 
krypton relative to deuterium were measured for the first 
time. The data show that the multiplicity ratio R'Ij is 
reduced at low v and high z. Different multiplicity ratios 
were observed for various hadrons. In contrast to the sim- 
ilarity between positive and negative pions, a significant 
difference in is found between positive and negative 
kaons and a larger difference between protons and an- 
tiprotons. The different results for various hadrons may 
reflect differences in the modification of quark and anti- 
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quark fragmentation functions |l2| and/or in the hadron 
nucleon interaction cross sections. 

The hadron multiphcity is observed to be enhanced at 
high in the nuclear medium, showing evidence of the 
Cronin effect in deep-inelastic scattering process. This 
effect is similar to the one observed in hadron nucleus 
scattering, with a rise of R'^j to values above unity for 
p2 > 1 GeV2. 

Additional measurements of differential hadron mul- 
tiplicities on both light and heavy nuclei with pion, 
kaon and proton identification are underway at HER- 
MES. Such measurements will also help to clarify the 
issues raised by the present data concerning the mass- 
number dependence of the hadronization process and of 
the Cronin effect for various identified hadrons. 
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